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Recently observed Aharonov-Bohm quantum interference of period h/2e in charge density wave 
rings strongly suggest that correlated density wave electron transport is a cooperative quantum phe- 
nomenon. The picture discussed here posits that quantum solitons nucleate and transport current 
above a Coulomb blockade threshold field. We propose a field-dependent tunneling matrix ele- 
ment and use the Schrodinger equation, viewed as an emergent classical equation as in Feynman's 
treatment of Josephson tunneling, to compute the evolving macrostate amplitudes, finding excellent 
quantitative agreement with voltage oscillations and current- voltage characteristics in NbSes- A 
proposed phase diagram shows the conditions favoring soliton nucleation versus classical depinning. 

PACS numbers: 71.45.Lr, 75.30.Fv, 03.75.Lm, 74.50.+r, 72.15.Nj, 73.23.Hk 



Cooperative quantum tunneling has emerged as an im- 
portant class of phenomena, whose manifestations in- 
clude Josephson tunneling [l| , macroscopic quantum tun- 
neling, and decay of the false vacuum [2]. The latter 
describes instability of a scalar field 0(r), where "vac- 
uum" refers to a minimum energy state. In this Letter, 
(/) represents the phase of a density wave. If (j){r) sits in 
a met ast able well ('false vacuum') it is unstable to decay 
by tunneling into a lower potential well within a small 
region, nucleating a bubble of 'true vacuum' bounded by 
solitons [§1. Herein we propose coherent Josephson-like 
tunneling of microscopic quantum solitons (single-chain 
solitons, delocalized in both longitudinal and transverse 
directions) of charge ±2e within a quantum fluid, i.e. not 
macroscopic quantum tunneling (MQT) of a massive ob- 
ject. 

The charge density wave (CDW) exhibits a charge 
modulation pix^t) = po{x,t) + pi cos[2/ci?x — (j){x^t)] 
along the chain direction, while the spin density 
wave (SDW) is equivalent to two out-of-phase CDWs for 
the spin- up and -down subbands Like a supercon- 
ductor, the density wave (DW) is a correlated electron 
(or electron-phonon) system capable of collective charge 
transport. The superconducting condensate is a charged 
superfluid represented by a complex order parameter. Its 
behavior can be described by the Schrodinger equation as 
an emergent classical equation for the condensate (Ch. 21 
of [5]). Unlike a superconductor, however, the order pa- 
rameter corresponding to the DW charge or spin mod- 
ulation does not couple directly to an electric field or 
vector potential. Nevertheless, gradients or kinks in DW 
phase carry charge that: 1) couple to an externally ap- 
plied electric field, and 2) generate their own electric field 



that leads to a Coulomb blockade effect. 

We propose that nucleated droplets of many ±2e 
charged kinks and antikinks behave as quantum fluids 
due to interchain interactions and quantum delocaliza- 
tion. We use the time-dependent Schrodinger equation 
to describe the coupled macrostates in a manner, as in 
Josephson tunneling, that is classically robust against 
decoherence below the transition temperature. In the 
proposed Josephson-like tunneling process, the quantum 
fluid flows over or through the barrier over a long time 
scale (up to ~ 1/is). Thermal excitations are frozen out 
by the Peierls gap since the condensate has one thermal 
degree of freedom within a phase-coherent domain . 

Aharonov-Bohm oscillations of period h/2e in the 
CDW magneto-conductance of NbSes crystals with 
columnar defects [7] and TaSs rings [8] suggest cooper- 
ative quantum behavior, in some cases over distances of 
85 pm and for T > 77 K. Moreover, the /i/2e, rather 
than h/2Ne period predicted Q for N parallel chains, 
supports the idea of coherent Josephson-like tunneling of 
microscopic entities of charge 2e within a quantum fluid, 
rather than MQT of a massive object. 

While the classical DW depinning field Ed is well un- 
derstood, less widely known is the existence of a Coulomb 
blockade threshold field Et (smaller than Ed) above 
which the system becomes quantum mechanically unsta- 
ble [10]. This threshold is readily determined for nu- 
cleation of charge soliton pairs in 1-D [iqI-[i2| or (in 3- 
D) soliton domain wall pairs, of charge j=Qo = ±2A/'epc 
where pc is the condensate fraction [13|. Just like the 
charged electrodes of a parallel-plate capacitor, these pro- 
duce an internal field, = Qo/eA. If an external field E 
is applied, the difference in electrostatic energies with. 
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^e'^{E ± E*)'^, and without the pair, ^eE'^, is positive 
when \E\ < Et = ^E*, yielding the Griiner relation [^]: 
cEt ~ epcN/A [14]. When added to the periodic pin- 
ning energy, this quadratic electrostatic energy ensures 
that the DW phase sits in the lowest potential energy 
well, or 'true vacuum' state when |£^| < Et- However, 
when E > Et, OT 6 = 2itE/E* > tt, the formerly 'true 
vacuum' becomes a met ast able state or 'false vacuum' 
(Figs. [Ik &b). 

Density waves usually have anisotropic relative dielec- 
tric responses, ey ^ e^, vs. the chain direction. Using 
rescaled coordinates, x' = x/ey, etc., a single-chain dis- 
location pair looks like a parallel plate capacitor that 
produces a field ^* = 2e/2eAch, where Ach is the cross- 
sectional area of a DW chain and e = eyeo Q- Thus, 
the Coulomb blockade threshold is comparable to that 
for domain wall pair creation, within a factor of ~ 1/2 . 
Many nucleated 27r dislocations, of charge ±2e each [l5|, 
can then form droplets with quantum fiuidic properties. 
The temperature-dependence of Et goes inversely with 
that of e(T) [14]. 

Our model relates the 'vacuum angle' (e.g. ref. [iq| 
and citing papers) to displacement charge Q between con- 
tacts by: = 27r{Q/Qo). The potential energy of the k^^ 
chain can then be written as |10|, ll2| : 



u [^fe] = 2uo [1 - cos(/)fe(x)] ^ue[0 - (l)k{x)Y 



(1) 



where the first term is the periodic DW pinning energy. 
The quadratic term is the electrostatic energy resulting 
from the net displacement charge or, equivalent ly, the 
applied field and internal fields created by kinks due to 
phase displacements. Fig. [TJi plots u vs. 6 when the en- 
ergy is minimized for (j) ~ 27rn (dropping the subscript) 
when ue << uq. The phases (j)k tunnel coherently into 
the next well via a matrix element T (Fig. [TJd) as each 
parabola, or branch (Fig. [1^), crosses the next at the in- 
stability points = 27r(n + 1/2). Here we propose an ide- 
alized time-correlated soliton tunneling model to simu- 
late DW dynamics. It includes a shunt resistance i?, rep- 
resenting normal, uncondensed electrons, in parallel with 
a capacitive tunnel junction representing soliton tunnel- 
ing (Fig.[Th), by analogy to time-correlated singe-electron 
tunneling (SET) [16]. 

Advancing the phase of all parallel chains by 27rn cre- 
ates multiple pairs of soliton domain walls that quickly 
reach the contacts. Similar to SET, the voltage is then 
proportional to net displacement charge: V = {Q — 
nQo)/C = {Qo/27rC)[0 — 27rn] when the phase has ad- 
vanced to (0) = 27rn between the contacts. More gen- 
erally, if the phase expectation value {(j)) among N par- 
allel chains advances by a fraction or non-integer multi- 
ple of 27r, the voltage is F = {Qo/27rC)[6 — {(j))] where 
C = eA/£. This leads to a total current: I = In -\- Idw^ 
where In = {Qo/27rRC)[0 — ((/))] is the normal current, 
and Idw — 
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FIG. 1: (color online) a. Potential energy vs. for (j) ^ 27rn. 
b. 1^(0) when = 2nE/E* > tt as the phases (/)/c(x) tun- 
nel coherently into the next well. c. Time- correlated soliton 
tunneling model. 



uj = 27tI/Qo and r = RC yields the following equation 
for the time evolution of 0: 



de 



1, 



uj--[e-m. 



(2) 



We compute {(j)) by solving the Schrodinger equation. 
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dt 



1,0 



(3) 



is the DW current. Defining 



to compute the original and emerging macrostate ampli- 
tudes ^|Jo{t) and '0i(t) (more generally ipn and V^n+i) for 
the system to be on branches and 1 (or n and n + 1, 
Fig. [T^), respectively. We interpret these amplitudes to 
represent classically robust order parameters, and the 
above equation is viewed as an emergent classical equa- 
tion following Feynman [5[ . The macrostates are coupled 
via a tunneling matrix element T with a Zener-like field 
dependence. (Another approach, employing probabilities 
rather than amplitudes, yields sharp sawtooth-shaped 
voltage oscillations and will be discussed elsewhere.) 

Our model represents the amplitudes ipo^i by: ?/^o,i = 
y^pojexp [z^o,i], where po,i = Nq^i/N is the fraction 
of parallel chains on the respective branch. Advanc- 
ing 0/c(x) by 27r within a given region, taking (pk from 
one branch to the next, is equivalent to creating a pair of 
microscopic 27r-solitons. Thus, the macrostate order pa- 
rameters V^o,i are coupled via coherent, Zener / Josephson 
tunnehng of delocalized quantum solitons [17|, with an 
enormous aggregate of N (up to ~ 10^) such processes 
occurring coherently. 

The driving force is the energy difference per unit 
length between potential minima at ^ 27rn and ~ 
27r(n + 1). When a = ue/uq << 1, this force is given 
by: F = 47r2X£; l9' where 0'^ = - 27r(n + ^). Fol- 
lowing Bardeen [li, [TH, T is estimated as: T{F) = 
-4FAexp[-Fo/F], where A"^ ^ A^/hvo + A;;,^ A^ is 
a mean free path length, A^^ is the microscopic soliton 
energy, vq is the phason velocity, and Fq ~ A'^/hvQ. This 
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expression is similar to the rate of Schwinger pair pro- 
duction in 1-D [20] . Since any negative energy difference 
(Figs. [TJi & b) within the 'bubble' along the x- direc- 
tion is balanced by the positive soliton pair energy at 
its boundaries, T couples states of equal energy, /7o = 
Ui = U. Thus, defining ?/^o,i = Xo,i(^) exp[— z/7t/fi], the 
Schrodinger equation (Eq.(3)) reduces to: ifidxo,i/dt = 
TXo,i. 

We define: t' t/r, f ujt/2it (cx /), q e/2iT,qo = 
Fo/2Ft = 6>o/27r, Ft = 2eET, and q'^ = = q-n- 

|, to simplify the computations. Finally, setting Xo(^) = 
co(t) and Xi(^) = ^ci(t), taking cq and ci to be real, yields 
the coupled equations: dci/dt' = [7^^ ^xp(— go/^o)] 
and dco/dt' = — [7^^^ ^xp(-g'o/^n)] Q.'n > 0? where 

7 = 327r'^UE)^r/h. These are integrated numerically, with 
initial values cq = 1 and ci = 0, yielding: {(/)) = 27r[n+p], 
where p = |ci|^. The transition from branch n to n + 1 
is considered complete, and n is incremented while p is 
reset back to zero, once p exceeds a cutoff close to one 
(e.g. 0.9995). When an applied current pulse is turned 
off, any remaining displacement charge discharges back 
through the shunt resistance and the system retains a 
memory of the previous macrostate amplitudes. The 
algorithm thus incorporates backward transitions from 
branch n to branch n — 1 when dQ/dt^ F, and q'^ are neg- 
ative. We find that no more than three training pulses 
are needed to converge to a 'fixed point' of one or two 
voltage oscillation patterns, which are averaged. 

Fig. [2^ compares the quantum theory with measured 
voltage oscillations [21| of NbSes for rectangular current 
pulses. The parameters used for the theoretical plots 
(solid lines) in Fig. [2^ are: 7=0.5, qo=0.7^ r=51 ns, 
V* = E^i=l.ll mV (where i= distance between con- 
tacts), and Rn= 99.6 ^. The measured threshold current 
of 6.93 /iA is taken to correspond to /= 0.6, the normal- 
ized onset threshold current consistent with the chosen 
values of 7 and qo, while the remaining normalized cur- 
rent pulse amplitudes / are scaled to the indicated am- 
plitudes in Fig. [2^. 

The theoretical voltage oscillation amplitudes agree 
remarkably well with experiment, especially consider- 
ing the simplicity of the model. Moreover, the model 
correctly reproduces the observed progression of non- 
sinusoidal shapes, ranging from rounded backward saw- 
tooth behavior for the 9.90-/iA current pulse to more 
symmetrical oscillations for higher current pulse ampli- 
tudes. The oscillation frequency is / = Idw/Qo^ so the 
number of oscillations per pulse, captured correctly by 
the model, increases with current. During each cycle, the 
CDW current, Icdw = I — V/Rn, gradually increases over 
a significant portion of a cycle as the voltage decreases 
from its maximum value. As seen in the bottom plot of 
Fig. [2^, this time scale can be up to ~ 1 /is, supporting 
the idea that the quantum fluid flows through the barrier 
for a relatively long time. 

The I — V and dV/dl curves are computed by aver- 
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FIG. 2: (color online) a. Theoretical (solid lines) vs. exper- 
imental (dashed lines ^]) voltage oscillations of an NbSes 
crystal at 52 K for current pulse amplitudes (bottom to top, 
offset by 0, 0.25, 0.5, and 0.75 V for clarity): 9.90 /iA, 
10.89 /iA, 11-49 /xA, and 11.88 /xA. b. Simulated DW cur- 
rent vs. field for 7= 1.5 and several go- Dotted lines (Ta- 
ble I): Bardeen's modified Zener function [22[. c. Simulated 
R — dV/dl vs. ///*, where /* = E^'i/Rn , where Rn is 
the normal resistance at zero bias, for several qo and 7= 1.5. 
d. Theoretical (solid lines, Table [ij) vs. experimental (dotted 
lines) dV/dl vs. / for an NbSes crystal. 



aging the voltage over several cycles, with results shown 
in Fig. [JJd & c. A range of behaviors are captured, with 
rounded Zener-like behavior (e.g. Refs. |22|, |23|), emerg- 
ing for large go ^ Eq/Et^ as contrasted with, when go 
is small, more Hnear I — V curves and dV/dl curves 
with negative dips or wings, as seen in NbSes crystals 
with fewer impurities \2A\. For small go, the 'measured' 
threshold Erm occurs near the Coulomb blockade thresh- 
old: Etui = Et = E^ /2. However, Etui becomes larger 
than E^ /2 as go increases. The dotted lines in Fig. [JJ) 
are obtained from a normalized Bardeen function [22|, 
iDw/r = m - E^tJ^M-KJEI where E^ = 
E/E* and /* = EH/Rn, while E'j.^ = Erm/E* and 
£^0m — Eom/E"^ are normalized 'measured' threshold and 
Zener activation fields, and F = 1.0 for all three plots. 
The remaining parameters used for the Bardeen function 
fits are shown in Table H 

The dotted lines in Fig. [2]i show, for an NbSes crystal, 
differential resistance, R = dV/dl^ normalized to normal 
resistance Rn = dV/dI\zerobias vs. I/lTm^ where / is 
total applied current and Irm = VTm/Rn is the mea- 
sured threshold current. The solid lines in Fig. [2li are 
simulation results using the parameters indicated in Ta- 
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TABLE I: Parameters used to generate the Bardeen function 
plots in Fig. 



Qo 






2 


0.847 


0.96 


4 


1.10 


2.55 


6 


1.40 


4.05 



bleini Fig. [211 shows excellent agreement between theory 
and the dV/dl measurements, showing rounded behavior 
below the upper and lower Peierls transitions. 



TABLE IL Parameters used to simulate the solid dV/dl 
curves in Fig. [2J1. 



Temperature 


iTm/r 


7 




20 K 


0.87 


2.7 


3.0 


25 K 


0.76 


3.2 


2.3 


35 K 


1.27 


2.8 


6.7 


70 K 


1.71 


0.41 


8.5 


120 K 


2.22 


0.275 


10.0 



Some CDW crystals exhibit more than one thresh- 
old within certain temperature ranges [H, [26| (Fig. [3^). 
The two major thresholds emerge naturally, provided 
the nucleated soliton conductance is sufficiently small 
for to be treated quasi-statically, i.e. 9 — iieE / ciEt^ 
where ei = e{E ^ Et)- We interpret the low- and high- 
field thresholds as due to soliton nucleation and clas- 
sical depinning, respectively. Fig. [Sb (left) illustrates 
the quantum {0 > tt) and classical {0 > Oc) instabili- 
ties, where Oc{(^) — (^~^ +7r/2 when a = ue/uq << 1. 
Fig. [Sb (right) plots u{(j)) when ^ = tt for several values 
of a. Figure [3fc shows the resulting 6 vs. ue/uq phase 
diagram, which illustrates the pinned state, for ^ < tt 
and ue/uq < 1, a region {tt < < Oc) in which soliton 
nucleation occurs, and a high field classical depinning re- 
gion {0 > Oc). 

The observed flat ac responses [27| and small phase 
displacements [28|] below threshold in NbSes and TaSs 
suggest ue/uq << 1 (solid arrow, Fig.[3t), where soliton 
nucleation dominates. For example, computed [12] phase 
displacements (0) below threshold compare favorably to 
the measured 2° phase displacement obtained from NMR 
experiments [28| on NbSea, provided ue/uq is taken to 



be 0.015 [12|. Using, ue/uq = 2iiEt/Ecu the 48 K 
blue bronze data [25| in Fig. [3^1 suggests a similar value 
of about 0.01. However, the soliton nucleation thresh- 
old field Et^ and consequently ue/uq^ increase with de- 
creasing temperature, whereas the classical threshold Ed 
shows a weak temperature dependence. We interpret the 
change in soliton nucleation threshold (which scales in- 
versely with e) as due to a reduction in e as the normal 
carrier concentration goes down with decreasing temper- 
ature. At 4 K the normal carriers become frozen out. 




ia' ia' io° 
Voltage (V) 



FIG. 3: (color online) a. Blue bronze I — V curves [25l l. 
b (left): u vs. 0, showing > n quantum instability and 
> Oc classical depinning. b (right): 1^(0) at ^ = tt 
for several a — ue /uq. c. Phase diagram showing pinned, 
soliton nucleation, & classically depinned states. Solid ar- 
row: ue/uq « 1 for which soliton pair creation dominates. 
Dotted arrow shows both soliton nucleation and classical de- 
pinning. Since ue oc 1/e, the path curves to the left (right) 
if e increases (decreases). Dashed arrows (a and c): classical 
depinning dominates. 



resulting in a relatively low e and sufficiently high ue/uq 
and Et for classical depinning to dominate (dashed ar- 
rows in Figs. [3^1 and[3t). 

A potential topic of interest is to study coupling of 
static and dynamic vector potentials to the relative phase 
6i — So between order parameters. This will be criti- 
cal, both for understanding the CDW ring Aharonov- 
Bohm experiments \^ and for interpreting ac response 
experiments [l^ , which show remarkable agreement with 
photon-assisted tunneling theory. Another variation of 
the model represents multiple DW domains, due to ran- 
dom pinning, as a network of many resistively shunted 
junctions of type shown in Fig. [It. 

Density wave transport is one of very few known cases 
of correlated transport of macroscopic numbers of elec- 
trons over long distances — the only known example of 
large-scale collective electron transport at biological tem- 
peratures (e.g. NbSa, with Tpdcris ^ 360 K [13] )• It 
is hoped that this Letter will revitalize this important 
branch of condensed matter physics, for which quantum 
principles have largely been ignored by most for the past 
thirty years. Additional areas of impact include improved 
understanding of other correlated electron systems, flux 
vortex nucleation, tunneling in quantum cosmology 29 1, 
and = TT instabilities in spontaneous CP violation [31 1. 
Finally, understanding of the quantum behavior of soli- 
tons could potentially lead to topologically robust forms 
of quantum information processing. 
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Dislocation pair nucleation. Density waves are often highly anisotropic, where the 
dielectric response, 6||, along the chain direction is much greater than that, e^, in the per- 
pendicular directions. In this case, the Coulomb blockade threshold for nucleation of a 27r 
soliton-antisoliton (dislocation- antidislocation) pair (Fig. S-la) is comparable to that for do- 
main wall pair creation, within a factor of ^ 1/2. This is seen by starting with the Maxwell 
equation: V • D = p, where (using the summation convention): Di = coCijEj. Here e^j is the 
relative dielectric tensor, which is diagonal with elements e\\ = Cxx^ = ^yy and e±z = ^zz 
if the axes i^j = x,^, and z are along principal crystallographic directions. Introducing 
rescaled variables, = x/e\\^ y' = y/e±y^ and = z/e±zi yields: • E = p/eo, resulting in 
the modified charge distribution, electrostatic potential, and field lines shown in Fig. S-lb, 
which assumes that e|| >> e±x,y In the rescaled coordinate system of Fig. S-lb, the disloca- 
tion pair looks like a parallel plate capacitor that produces an internal field £"* = 2e/2eAch^ 
where Ach is the cross-sectional area of a DW chain and e = e\\eo. Fig. S-lc illustrates the 
aggregation of many 27r dislocations of charge 2e into soliton domain walls, each of which 
may behave as a quantum fiuid and between which the bubble of 'true vacuum' grows. 

Impurity and temperature dependences of quantum and classical threshold 
fields. The Coulomb blockade threshold of a density wave scales inversely with low fre- 
quency dielectric response e which, in turn, is strongly infiuenced by pinning due to impurity 
concentration n^, normal carrier concentration, and temperature T. The Griiner relation be- 
tween dielectric response and threshold, for the simplest classical model, is given by 
eEci = AnerichPc-) where rich = N/A and Ed is the classical depinning field. The Coulomb 
blockade threshold Et = Qo/eA leads to a similar relationship: cEt = erichPc^ with the 
result, Et = E^/^'n^ which scales as for a density wave weakly pinned by impurities. 
Screening by normal carriers, both thermally activated and due to incomplete Peierls gaps 
in materials such as NbScs, further enhance e and reduce the ratio Et/E^i. For a fixed n^, 
the temperature dependence of e{T) leads (inversely) to the observed temperature depen- 
dence of Et- This leads to rather complex temperature dependence in materials such as 
NbSea, for example, which undergoes two separate Peierls transitions. 

The possibility of additional I — V infiection points at intermediate fields in the soliton 
nucleation region below the classical depinning field is suggested by some experiments on 
TaSs Rl? and might occur if the first nucleated solitons were impeded from rapidly reaching 
the contacts, perhaps by disorder. A second bubble bounded by solitons could then nucleate 
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within the first bubble. The additional electrostatic energy would be: ^eA£[{E—2E*'Y — {E— 
E"*')^] = eAlE'^'WE'"' — E]^, which is positive unless E exceeds a higher Coulomb blockade 
threshold, Et?^ = = S^Eti- (Note that, in general, 63 = e{E = Ets) ^ eie{E = En).) 
In a similar fashion, a series of higher Coulomb blockade thresholds Etu could exist, provided 
disorder or other mechanisms prevented solitons from annihilating or reaching the contacts 
quickly. Furthermore, 47r, 67r, etc. soliton domain walls might nucleate, presumably with 
much lower probability, causing additional inflection points at ^ = 27rn. A generalization 
of the Griiner relation would then have the form: e^ETn = '^^'^chPc-> where n is an integer. 
Such possibilities can potentially be explored in crystals with both weak and strong pinning 
centers, perhaps through ion irradiation of initially weakly pinned samples. 

Discussion of relevant energy scales. The microscopic quantum soliton energy 
gap (per dressed electron pair) A^^ can be estimated from the measured Zener fleld £"0, 
typically ^ 10 V/m. Using a phason velocity v^) = {m/Mjp)^/'^VF ^ 3 x 10^ m/s, where 
Mf/tti is the Frolich mass ratio , and Eo /^l/Tiv^e, yields - 4 x 10"^ eV. Al- 
though this microscopic energy is extremely small, A^^ <C fcT, both the initial and emerging 
macrostates become macroscopically occupied with substantial condensation energies due 
the large number (^ 10^) of interacting parallel chains. Thus the condensed quantum soli- 
tons and antisolitons in the emerging macrostate are effectively trapped in quantum fluids, 
preventing thermal excitations except across the Peierls gap, Ap > /cT, typically several 
lO's of meV. 

Two analogies are helpful here. One is coherent Josephson coupling between two super- 
conducting order parameters that are superimposed in space, as discussed by Leggett [4]. 
In this case, it is immaterial whether the macrostates are separated by energies per electron 
pair that are small compared to /cT, since thermal excitations only occur across either of 
the two BCS energy gaps. Another analogy is that of quantum creep of superfluid helium 
out of a container. If the height of the rim is, for example, d ^1 cm above the liquid surface 
inside the container, then the gravitational energy barrier for each helium atom, mgd is only 
^ 4 X 10~^ eV, which is small compared to kT even for temperatures as low as 1 mK. Nev- 
ertheless, the superfluid helium atoms remain trapped in the quantum fluid and prevented 
by the condensation energy from thermally hopping over the barrier as individual particles. 



3 




FIG S- 1: (a). COMSOL Multiphysics (COMSOL, Inc.) simulation of the electrostatic potential 
(red = positive, blue = negative) and electric field lines for an electric dipole consisting of dis- 
locations represented as positive and negative rectangular charge distributions, (b). COMSOL 
simulation for the same system, but with highly anisotropic dielectric constants, in the rescaled 
coordinate system. The apparent distance between the charges is greatly reduced along the x' 
direction, yielding an arrangement resembling a parallel plate capacitor, (c). In the original coor- 
dinate system, aggregation of many dislocations into 'domain walls', between which the bubble of 
'true vacuum' grows as the domain walls are driven toward the contacts by the externally applied 
electric field. 



* U hmiller@uh.edu 
[1] G. Griiner, Reviews of Modern Physics 60, 1129 (1988). 

[2] G. Griiner, Density Waves in S'o/zrf5( Addison- Wesley, Reading, Massachusetts, 1994), p. 259. 
[3] S. V. Zaitsev-Zotov, G. Remenyi, and P. Monceau, Physical Review Letters 78, 1098 (1997). 
[4] A. J. Leggett, Progress of Theoretical Physics 36, 901 (1966). 



4 



